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ABSTRACT 

We report on a 39 ks observation of the 7.7-s low-mass X-ray binary pulsar 4U 1626-67 with the High Energy 
Transmission Grating Spectrometer (HETGS) on the Chandra X-Ray Observatory. This ultracompact system 
consists of a disk-accreting magnetic neutron star and a very low mass, hydrogen-depleted companion in a 42-min 
binary. We have resolved the previously reported Ne/O emission line complex near 1 keV into Doppler pairs 
of broadened («2500 km s -1 FWHM) lines from highly ionized Ne and O. In most cases, the blue and red line 
components are of comparable strength, with blueshifts of 1550-2610 km s" 1 and redshifts of 770-1900 km s -1 . 
The lines appear to originate in hot (« 10 6 K), dense material just below the X-ray-heated skin of the outer 
Keplerian accretion disk, or else possibly in a disk wind driven from the pulsar's magnetopause. The observed 
photoelectric absorption edges of Ne and O appear nearly an order of magnitude stronger than expected from 
interstellar material and are likely formed in cool, metal-rich material local to the source. Based on the inferred 
local abundance ratios, we argue that the mass donor in this binary is probably the 0.02 Mq chemically fractionated 
core of a C-O-Ne or O-Ne-Mg white dwarf which has previously crystallized. 

Subject headings: accretion, accretion disks — binaries: close — stars: individual (4U 1626-67) — stars: 
neutron — techniques: spectroscopic — X-rays: stars 



1. INTRODUCTION 

The ultracompact low-mass X-ray binary 4U 1626-67 con- 
sists of a 7.7 s pulsar accreting from an extremely low-mass 
(< O.IMq) companion with a binary separation of only sd lt- 
sec. The very short binary period (P or b=42 min; Middleditch 
et al. 1981; Chakrabarty 1998) indicates that the companion 
must be extremely hydrogen-poor (Paczynski & Sienkiewicz 
1981; Nelson, Rappaport, & Joss 1986). Based on X-ray tim- 
ing limits, the Roche-lobe-filling companion must be either a 
0.02 Mq white dwarf or a 0.08 Mq hydrogen-depleted, partially 
degenerate star (Levine et al. 1988; Verbunt, Wijers, & Burm 
1990; Chakrabarty 1998). These solutions correspond to rel- 
atively face-on binary inclinations of 33° and 8°, respectively. 
For a randomly oriented ensemble of binaries, the a priori prob- 
abilities of finding such low inclinations are only 16% and 1%, 
respectively 2 . 

The pulsar's spin frequency evolves on short (\v jv\ ~ 6000 
yr) time scales due to accretion torques. After spinning up 
steadily for 13 years, the accretion torque abruptly changed 
sign in 1990, leaving the pulsar in a prolonged spin-down state 
(Chakrabarty et al. 1997). The surface dipole magnetic field 
strength of 3 x 10 12 G, inferred from the detection of cyclotron 
absorption lines (Orlandini et al. 1998), suggests that the pulsar 
is near spin equilibrium. The ultraviolet spectrum of the source 
is consistent with that of an accretion disk whose inner edge is 
truncated by the pulsar's magnetosphere near the corotation ra- 
dius r co = (GM x P 2 pin /47r 2 ) I//3 = 6.5 x 10 8 cm, and whose outer 

edge is truncated at the binary's tidal radius r t w2x 10 10 cm 
(Wang & Chakrabarty 2001). 

The pulsar's X-ray continuum spectrum is evidently corre- 



lated with its torque state. During the 1977-1990 spin-up 
epoch, the pulsar's pulse-phase-averaged X-ray spectrum was 
roughly described by the sum of three separate power-law 
components with photon indices 3 71 w 1.4-1.65 (0.7-10 keV), 
72 « 0.5 (10-20 keV), and 73 w 5 (20-60 keV) (Pravdo et al. 
1979; Maurer et al. 1982; Eisner et al. 1983; Kii et al. 1986; 
Mavromatakis 1994). However, observations acquired after the 
1990 torque reversal measured a significantly flatter soft X-ray 
spectrum with 71 w 0.6-0.8 in the 1-10 keV range (Angelini 
et al. 1995; Owens et al. 1997; Vaughan & Kitamoto 1997). 
Since the discovery of this X-ray source over two decades ago, 
its overall X-ray flux has declined steadily (Chakrabarty et al. 
1997 and references therein; Owens et al. 1997), but the optical 
flux from the accretion disk has remained essentially constant 
over the same time interval (Chakrabarty 1998). 

X-ray spectroscopy with ASCA and BeppoSAX also de- 
tected a strong, unusual complex of emission lines near 1 keV; 
these were identified as the Lya lines of hydrogenic and He- 
like Ne and O (Angelini et al. 1995; Owens et al. 1997), with 
possible contributions from Fe-L shell emission. Angelini et al. 
(1995) also detected a faint feature near 1.41 keV, which sug- 
gests the presence of photoionized plasma if the feature is inter- 
preted as the recombination continuum of hydrogenic Ne. Inde- 
pendent of whether the line complex arises from a photoionized 
or a collisionally ionized plasma, however, its strength indicates 
an overabundance of both Ne and O relative to solar values (An- 
gelini et al. 1995). 

In this paper, we present high-resolution spectra of 4U 
1626-67 obtained with the High Energy Transmission Grating 
Spectrometer (HETGS; Canizares et al. 2001, in preparation) 
on board the Chandra X-ray Observatory. 



1 Also Department of Physics, Massachusetts Institute of Technology, Cambridge, MA 02139. 

2 A more massive 0.6 Mq helium-burning star is also allowed by the X-ray timing limits (with i =1.3° and an a priori probability of 0.03%), but this solution requires 
an implausible (36 kpc) distance to be consistent with the observed X-ray flux. 

3 We define photon index 7 such that the photon flux dN/dE oc E^ 1 . 
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2. CHANDRA OBSERVATIONS AND DATA REDUCTION 

4U 1626-67 was observed with the HETGS on 2000 Septem- 
ber 16 (05:47:21 UT) continuously for 39 ks. The HETGS car- 
ries two different types of transmission gratings: the Medium 
Energy Gratings (MEG) and the High Energy Gratings (HEG). 
These allow for high-resolution X-ray spectroscopy in the 1-35 
A (0.35-12.4 keV) range, with a spectral resolution of A/AA ~ 
1400 at 12 A (1 keV) and A/AA « 180 at 1.8 A (6.9 keV). 
The dispersed spectra were recorded with an array of 6 charged 
coupled devices (CCDs) in the focal plane which are part of the 
Advanced CCD Imaging Spectrometer (ACIS; Garmire et al. 
2001, in preparation). We defer to the Chandra X-ray Center 
(CXC) documents for a more detailed description of the spec- 
troscopic instruments (see http://chandra.harvard.edu). 

We recorded a total of 62373 events in the 1st order MEG 
spectrum and 33744 events in the HEG spectrum after applying 
standard event grade selection cuts. For some of the analysis, 
we combined the MEG and HEG spectra where their passbands 
overlapped. In these cases, we used a count-weighted average 
of the resolution of both gratings, which is between 0.01 and 
0.02 A. The CXC provided aspect-corrected "level 1" event lists 
via standard pipeline processing. These data were reprocessed 
using the latest available data processing input products. We 
also removed all events that were attributed to bad pixels, bad 
columns, or flaring background events that were not handled by 
the CXC standard processing. 

The determination of the zeroth order source position is cru- 
cial for the calibration of the wavelength scales, since it defines 
the zero point for wavelength in the dispersed spectra. We can 
determine the source position on the detector to an accuracy of 
0.2 pixels, which corresponds to a wavelength zero point accu- 
racy of 0.002 A in MEG and 0.001 A in the HEG 1st order. The 
current status of the overall wavelength calibration is of order 
~ 0.05%, leading to a worst-case uncertainty of 0.004 A in the 
1st order MEG at O VII (21.6 A), and 0.006 A in the 1st order 
HEG at Ne X (12.1 A). We processed the event lists into "level 
1 .5" grating event lists using available CXC software. For the 
spectral analysis we used a combination of custom software, 
FTOOLS, and the CIAO 2.1 package from CXC. We computed 
aspect-corrected exposure maps for each spectrum, allowing us 
to correct for effects from the effective area of the CCD spec- 
trometer. The overall flux calibration throughout the CCD array 
is of order 10%, although it can be as poor as 20% at lower en- 
ergies in some areas. 

The zeroth order image was affected by heavy "pileup", 
where the event rate was so high that two or more events were 
detected in the CCD during a 3.2 s frame exposure. The time 
between frames was dt = 3.24104 s, which accounts for the 
time to transfer charge to the frame store. Pileup distorts the 
CCD spectrum because detected events overlap so their de- 
posited charges are collected into single, apparently more en- 
ergetic events. More significant, however, is that many events 
are lost on-board as the grades of the piled events overlap those 
of rejected background. In this case ~ 90% of the events are 
lost. We ignore the zeroth order events in all subsequent analy- 
ses. 

MEG and HEG events were selected for timing analysis on 
the basis of spatial, pulse height, and grade selection (as de- 
scribed in the last section). Events from first and second or- 
ders were included, giving 102,835 events. For the same pulse 
height and grade selection, a background region 12.6" away 



gave only 386 events in the same spatial window, so back- 
ground is ignored. Event times were corrected for the motion 
of the spacecraft and then to the solar system barycenter. A ran- 
dom quantity uniformly distributed between and dt was added 
to event times in order to avoid aliasing with the sample time. A 
significant peak was found in the FFT at a period of 7.672 s and 
a x 2 period folding method gave a best period of 7.6717(7) s, 
consistent with an extrapolation of the source's long-term spin- 
down behavior (Chakrabarty et al. 1997). The pulse appears 
sinusoidal with an amplitude of 4.0 ± 0.9%, averaged over the 
entire bandpass. This sinusoidal appearance may as well be a 
consequence of the fact that the pulse period is very close to the 
time resolution of the detector. It has been shown in previous 
observations that the shape and amplitude of the intrinsic pulse 
chnage with energy (Pravdo et al. 1979, Kii et al. 1986, Owens 
et al. 1997). In this respect we observe that the modulation 
of the events with E > 2 keV is 5.3 ± 1.3%, somewhat higher 
than the modulation found in the lower energy band (E < 2 
keV): 2.7 ±1 .2%. Above 2 keV, the hard power law dominates 
the spectrum (see section 4.2) while below 1 keV, the power 
law and blackbody components have about equal contributions 
(see Fig. 4), so this marginally significant change of the pulse 
fraction could be obtained if the blackbody component is totally 
unpulsed. 

3. SPECTRAL ANALYSIS 

The preparation of the spectra for analysis was done in sev- 
eral steps. In order to perform the basic spectral fits, we first 
subtracted the background from each spectrum and fitted the 
HEG and MEG spectra simultaneously using the XSPEC vl 1.0 
package (Arnaud 1996). We also co-added all background- 
subtracted spectra and exposure maps and divided the spectra 
by these exposure maps in order to screen the fit against local 
deviations and adjust the fit. The exposure map consists of the 
aspect-corrected effective area of the instrument over the total 
exposure for each dispersion arm. The source spectra were ex- 
tracted within a cross-dispersion width of 2x 10~ 3 degrees en- 
compassing 98% of the dispersed flux on the CCD. We then 
selected the first-order spectra from the HEG and MEG. The 
background spectra were computed by extracting spectra above 
and below the dispersed flux, again with the same total width in 
cross-dispersion. The background usually appeared extremely 
weak, although we observed some enhanced flux in the MEG 
-1st order range between 16 and 20 A, at a level about a fac- 
tor of 2.5 above the average. Figure 1 shows the background- 
subtracted count spectrum after co-adding the MEG and HEG 
spectra. The lower curve shows the total contributing back- 
ground above 10 A. 

3.1. Photoelectric Absorption 

Photoelectric absorption by neutral metals along our line of 
sight through the interstellar medium attenuates the low-energy 
X-ray continuum. With HETGS, we are able to resolve the indi- 
vidual photoelectric absorption edges of Ne-K (14.3 A), Fe-Lni 
(17.5 A), and O-K (23.3 A), allowing separate determinations 
of the column density for each of these elements (see Table 1). 
Figure 2 shows the Ne K and Fe L-edges. We note that the 
uncertainties in these edges are still large, since the derived ab- 
solute edge depths are of the order of the statistical uncertainties 
of the data bins, and their determination therefore also depends 
on the proper placement of the continuum (see below). 
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FIG. 1. — The Chtmdra/HETGS count spectrum of 4U 1626-67. The upper 
curve is the total background-subtracted spectrum, and the lower curve is the 
background above 10 A. 

For Ne and O, we used the atomic cross sections of Verner 
et al. (1993) but for Fe-L we instead used more recent labo- 
ratory atomic constants (Kortright & Kim 2000; see also dis- 
cussion in Schulz et al. 2001). Also given in Table 1 is the 
column density of neutral hydrogen (Nh) implied by each ele- 
ment's absorption edge, assuming the abundances of Morrison 
& McCammon (1983) for the material along the line of sight. 
It is interesting to note that for Ne and O, these implied val- 
ues of N H are higher than the value of N H = (6.2 ± 0.7) x 10 20 
cm" 2 measured directly from Lycv measurements in the ultra- 
violet (Wang & Chakrabarty 2001). The abundance ratios also 
differ substantially from solar values. 

Although the threshold for the C-K edge (43.7 A) lies beyond 
our passband, the short wavelength "tail" of this edge does af- 
fect the continuum we observe. One way to obtain a better 
representation of the continuum in the wavelength range above 
18 A is to assume a C/H number ratio (5-1 8) x solar. For our 
continuum representations in the figures below, we have conser- 
vatively assumed a carbon overabundance of (C/H)/(C/H)q=5. 
Given that Nh ^3 x 10 21 cm" 2 , we expect no sensitivity to ab- 
sorption edges of Mg and higher-Z elements. However, it is 
interesting to note the presence of a marginal edge-like feature 
near the position of the Mg K-edge at 9.5 A. In order to explain 
this feature by Mg absorption, we would require an order of 
magnitude overabundance of Mg relative to the solar value. 



Table 1 

Photoelectric Absorption Towards 4U 1 626-67 



A 




Optical 


Total N z * 


Implied Nh 


Local Nz" 


(A) 


Edge 


Depth 


(10 17 cm" 2 ) 


(10 21 cm" 2 ) 


(10" cm" 2 ) 


9.5" 


Mg K (?) 


0.094 ± 0.052" 


3.8±2.l" 


9.4 ±5.2" 


3.5±2.l" 


14.3 


NeK 


0.19 ±0.09 


5.2±2.5 


3.7 ± 1.8 


4.4 ±2.5 


17.5 


Fe Lin 


0.25 ±0.11 


0.35 ±0.16 


1.1 ±0.5 


0.15±0.15 


23.3 


OK 


1.40±0.67 


25 ±12 


3.3 ± 1.6 


20 ±12 


43.7 


cr 




10-36" 


2.2-8. l e 


7-33 e 



a Total column density for specified element. 

b Nu implied by total Nz assuming standard ISM abundances (Morrison & McCam- 
mon 1983). For comparison, the directly measured Nh from Lya is (0.62 ± 0.07) x 
10 21 cnT 2 . 

"Column density local to source for specified element, after correcting for inter- 
stellar absorption using Nh = 6.2 x 10 20 cm -2 and stanndard ISM abundances. 

"Feature not reliably identified in our data. 

"Edge lies outside our bandpass. Strength is inferred indirectly from continuum 
strength redward of O edge. 



3.2. Continuum Model and Fluxes 

To fit the continuum spectrum, we rebinned the data into 
broad (0.08 A) bins, resulting in about 300 wavelength chan- 
nels. To avoid skewing of the fit due to obvious bright lines, 
we always included several Gaussian line components in the 
model fit. This task appeared to be quite difficult, since we 
observe quite a number of lines (see below) in the spectrum, 
which appear very broad and a substantial uncertainty remains 
due to the fact that we cannot identify possible contributions 
from many more weak broad lines. Spectral models with a sin- 
gle absorbed continuum component did not fit the data. How- 
ever, absorbed two-component fits were successful using either 
a pair of power-law components or a sum of power-law and 
blackbody components. The best-fit parameters are given in 
Table 2. Although the HETGS data does not clearly prefer one 
model over the other, the PL+BB model is consistent with pre- 
vious low-resolution broadband measurements (Angelini et al. 
1995, Owens et al. 1997), and we will adopt this model with 
the photon index fixed at 0.65 for the remainder of this paper. 
The spectrum versus wavelength is shown in Figure 3. A plot of 
the photon spectrum versus energy is shown in Figure 4, with 
the individual continuum components clearly indicated. 

The total measured X-ray flux is 2.09 xlO" 10 
erg cm" 2 s" 1 (0.5-10 keV). This corresponds to an (unab- 
sorbed) luminosity of 2.6 xlO 34 erg s" 1 d\ , which is similar 
to that reported by Owens et al. (1997). This would seem to 
break the trend of steady dimming over the entire source history 
(Chakrabarty et al. 1997). The blackbody temperature is also 
consistent with that measured by Owens et al. (1997). 

3.3. Emission Lin es 

In order to study discrete line features, we subtracted the con- 
tinuum model from the measured photon spectrum. This resid- 
ual spectrum is shown at the bottom of Figure 3. It shows a 
variety of emission and absorption features. We have chosen 
not to analyze the wavelength range below 8 A, since the pos- 
sible features are weak and are not easily disentangled from 
calibration uncertainties. The residual spectrum contains sev- 
eral strong emission line features at higher wavelength, which 
we identify exclusively from hydrogenic and He-like Ne and O 
ions. Fitting single Gaussians to the lines yielded line widths 
ranging from 0.25 A at Ne X to 0.48 A at O VII, corresponding 
to Doppler velocities of 5340 and 7460 km s" 1 , respectively. 
The centroid of the Ne IX lines is closest to the expected posi- 
tion of the intercombination line, while the O VII line lies be- 
tween the expected positions for the intercombination and res- 
onance lines. 

Single Gaussian line profiles did not provide an acceptable 
fit to the observed line shapes. Figure 5 shows the profiles of 
the neon and oxygen lines, using spectral binnings of 0.015 A 
and 0.06 A. The H-like lines clearly show a double line pro- 
file. We therefore fit all observed line profiles with Gaussian 
line doublets. The resulting line parameters are given in Table 
3. For the He-like triplet lines, we assumed that the individual 
components of each triplet were subject to the same velocity 
shift, and we also enforced an upper limit on the line flux using 
the 3 <j error level of the continuum. We also can further set a 
lower limit for these fluxes from calculations for photoionized 
plasmas by Porquet and Dubau (2000); specifically, the lines of 
both Ne IX and O VII should obey G < 4.9 (see below). 

For each line pair, we denote the higher and lower wave- 
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FIG. 2. — The spectrum of 4U 1626-67 near the Ne AT-edge (left) and the Fe L-edges (right). The data have been smoothed for illustrative purposes. The dotted 
lines show the expected contribution of interstellar photoelectric absoprtion based on an ultraviolet measurement of Nu ■ The solid lines show the best-fit model for 
the much stronger absorption edges in the Chandra data. 



length components as "red" and "blue", respectively, under the 
assumption that the wavelength offsets are caused by Doppler 
redshifts and blueshifts. Table 3 includes the measured wave- 
lengths and shifts for each of the observed lines. We observe 
blueshifts of 1500-2600 km s" 1 and redshifts of 800-1900 km 
s" 1 . The FWHM line widths for both the red and the blue com- 
ponents are roughly 2500 km s -1 for all lines. The equivalent 
widths for the integrated red+blue line pairs are 22.5 eV for 
Ne X, 17.8 eV for Ne IX, 24.6 eV for O VIII, and 44.4 eV for 
O VII. 

3.4. Plasma Diagnostics 

The relative strengths and shapes of the line components in 
a He-like triplet are a sensitive diagnostic of plasma temper- 
ature and density (Gabriel & Jordan 1973; Porquet & Dubau 
2000). We denote the resonance, intercombination, and forbid- 
den components as r, i, and /, respectively. We also denote 
hydrogenic Lyman a lines as h. Plasma temperature may be 
estimated from the ratio G = (f+ i)/r. At high (> 10 6 K) tem- 
peratures, there is an additional dependence on the ratio of hy- 
drogenic and He-like ions, H = h/r. Plasma density may be 
inferred from the ratio R = f/i. We compute these ratios for 
each of our red and blue line components separately, and then 
compare them with the calculations of Porquet & Dubau (2000) 



to infer plasma temperature and density. The results are sum- 
marized in Table 4. 

We infer temperatures of (1.6-1. 8)x 10 6 K from the Ne IX 
lines and (0.9-1 .0) x 10 6 K from the O VII lines. For our density 
estimates, the /lvalue calculations of Porquet & Dubau were 
in fact made for temperatures close to these, although our R- 
values are low enough to be fairly insensitive to temperature in 
any case. The density of the Ne IX regions is about an order of 
magnitude higher than in the O VII regions. 

Our observation of r < i for the He-like lines would normally 
be interpreted as indicative of a purely photoionized plasma, 
in contrast to a hybrid plasma which experienced collisional 
ionization as well; however, this seems surprising for the high 
densities we measure. The calculations of Porquet and Dubau 
(2000) do not include scattering effects, which may be impor- 
tant in high-density regimes. From the fact that G > 1 is well 
established for all of our line components, we can deduce that 
the emission arises from a recombining plasma (A. Pradhan 
2001, personal communication; Liedahl et al. 2001). Given 
that we observe moderately high densities in a hot recombining 
plasma, it seems likely that collisional ionization is indeed still 
important. The weak resonance line in the He-like triplets may 
therefore be an opacity effect due to enhanced line scattering. 



Table 2 
Continuum Spectral Fits 





Ah 


Power law 


Power law 


Blackbody 




Model 


(10 21 cm" 2 ) 


At a 


71 


A 2 a 


72 


RL/D z 10kpc kT (eV) 


xl 


PL+BB 
PL+BB 


0.26 ±0.06 
0.62 (fixed) 


8.93±0.14 
7.42 ±0.14 


0.84 ±0.07 
0.65 (fixed) 






55.7 ±25.1 333 ±50 
173.0 ± 16.0 336 ±42 


1.03 
1.05 


PL+PL 
PL+PL 


0.87 ±0.09 
0.62 (fixed) 


9.91 ±0.41 
6.77 ±0.84 


0.84 ±0.14 
0.65 (fixed) 


3.99 ±0.64 
5.87 ±0.92 


2.51 ±0.63 
1.74±0.30 




1.02 
1.06 



"Normalization of PL component at 1 keV in units of 10 3 erg cm 2 s 1 keV 1 . 
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FIG. 3. — The HETGS photon spectrum of 4U 1626-67 as a function of wavelength, along with the best-fit power law (dashed line) and blackbody (dotted line) 
components. The residual spectrum in the bottom panel identifies various broad emission lines as well as ionization edges. 
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FIG. 4. — The broad-band photon spectrum of 4U 1626-67 as a function 
of energy, showing the relative contributions of the blackbody (dotted line) 
and power law (dashed line) components. The power law component clearly 
dominates at high energy. 

This introduces some additional uncertainty to our measured 
G-value and the inferred temperature. 

4. DISCUSSION 

4. 1 . Origin of the Double-Peaked Emission Lines 

We have resolved the previously detected Ne/O emission line 
complex into Doppler pairs of broadened lines from hydrogenic 
and He-like Ne and O. For most of the line pairs, the red and 
blue components have comparable strength (the one exception 
being the Ne X line). X-ray line complexes with well-separated 
red/blue line pairs have been previously detected from relativis- 
tic jets in the so-called Galactic microquasars SS 433 (Kotani et 
al. 1996; Marshall, Canizares, & Schulz 2001) and IE 1740.1- 
2947 (Cui et al. 2000), but the Doppler shifts in those cases are 
relativistic, quite unlike what we observe in 4U 1626-67. In- 
stead, the velocity width of the lines and their projected Doppler 
velocity shifts are of the same order of magnitude as the Keple- 
rian velocities in the accretion disk (970-5400 km s" 1 ). Thus, 
the most natural explanation for these lines is that they arise in 
or near the Keplerian disk flow around the neutron star. Sim- 
ilar double line profiles are frequently observed in optical/UV 
emission lines from accretion disk flows in cataclysmic vari- 
ables (see, e.g., Frank, King, & Raine 1992), but our observa- 



tion would be the first X-ray detection of such features. (Single- 
peaked lines with comparable broadening and He-like line ra- 
tios have been observed in the edge-on LMXB EXO 0748-676; 
Cottam et al. 2001). We note that if this Keplerian explana- 
tion is correct, then the companion must be a 0.02 M Q white 
dwarf. The other alternative (a hydrogen-depleted, partially de- 
generate 0.08 M Q star) requires a very small (< 8°) inclination. 
This gives a maximum projected velocity at the inner disk edge 
of only < 750 km s" 1 , inconsistent with our observed Doppler 
shifts. 

Alternatively, if the companion really is a more massive star, 
then we are viewing the disk nearly face-on and the Doppler 
line pairs are more difficult to understand. One possible ex- 
planation is that they arise in a bipolar outflow. Several au- 
thors have proposed that a magnetohydrodynamic wind might 
be driven from the magnetopause (near the inner disk edge) of 
an accreting neutron star during spin down (Anzer & Borner 
1980; Arons et al. 1984; Lovelace, Romanova, & Bisnovatyi- 
Kogan 1995). Such an outflow would presumably have roughly 
the Keplerian velocity at the inner disk edge but would be ori- 
ented at a substantial angle to the disk plane. An advantage of 
this picture is that there is already indirect evidence that some 
outflow from the disk must be present, since the X-ray flux 
from the source (which traces M onto the neutron star) has been 
steadily dropping for years while the optical flux (which traces 
M through the disk) has remained unchanged. We can estimate 
the wind mass loss rate as -M w ind ~ An e m p r 2 V ~ 10" 10 M Q 
yr -1 , where we have used the inner disk radius and Keplerian 
velocity and have taken A = 16 for oxygen. This is of the same 
order as the observed mass accretion rate onto the neutron star. 
However, a serious drawback of this explanation is that it fails 
to account for the comparable strengths of the red and blue line 
components, since we would expect the redshifted lines to be 
at least partially blocked from our line of sight by the accretion 
disk itself. Consequently, we do not regard this explanation as 
likely. 

The He-like triplets of O and Ne are dominated by the in- 
tercombination component, indicating that the lines arise in a 
relatively dense, photoionized plasma (see, e.g., Liedahl et al. 
2001). The inferred plasma temperatures (see Table 4) are con- 
siderably higher than the expected surface temperature range 
of the optically thick accretion disk. However, detailed calcu- 



Table 3 



Emission Line Parameters 











Blueshifted lines 






Redshifted lines 








FWHM 




A\ c 


V 






V 


(A) 


Ion 


(kms- 1 ) 


Flux b 


(A) 


(kms" 1 ) 


Flux b 


(A) 


(kms" 1 ) 


12.130 


NeX 


2860 ±330 


8.15±0.93 


-0.090 ±0.014 


-2220 ±350 


15.04±1.65 


0.050 ±0.009 


1240 ±220 


13.440 


Ne IX (r) 


2320 ±300 


1.68±0.64 


-0.070 ±0.024 


-1550 ±530 


1.91±0.60 


0.035 ±0.013 


770 ±290 


13.550 


Ne IX (0 


2300 ±290 


5.19±1.79 


-0.070 ±0.024 


-1550 ±530 


5.41±2.00 


0.035 ±0.013 


770 ±290 


13.700 


Ne IX (/) 


2270 ±280 


1.69±0.64 


-0.070 ±0.024 


-1550 ±530 


1.92±0.62 


0.035 ±0.013 


770 ±290 


18.970 


OVIII 


3330 ±380 


14.04±2.52 


-0.1 10 ±0.028 


-1740 ±440 


17.82±0.57 


0.120 ±0.030 


1900 ±480 


21.600 


O VII (r) 


2430 ±310 


4.14±1.84 


-0.190 ±0.073 


-2610 ±1000 


5.26±2.48 


0.100 ±0.041 


1370 ±560 


21.800 


O VII (0 


2410 ±300 


9.02±3.64 


-0.190 ±0.073 


-2610 ±1000 


9.59±3.83 


0.100 ±0.041 


1370 ±560 


22.100 


OVII(/) 


2380 ±300 


2.07±1.03 


-0.190 ±0.073 


-2610 ±1000 


1.88±0.98 


0.100 ±0.041 


1370 ±560 



"Calculated values from Mewe et al. (1985). 
b l(r 5 photons cm^s -1 
c Measured shift from (a) in MEG. 
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FIG. 5. — Double-peaked, broad emission lines from 4U 1626—67. The upper panels show the hydrogenic Ne X (upper left) and O VIII (upper right) lines, and 
the lower panels show the He-like triplets of Ne IX (lower left) and O VII (lower right). The blue curves denote the model for the blue-shifted components, the red 
curves the model for the red-shifted components, and the green curves the model for the overall line shape. In addition to the Doppler shifts, the individual line 
components are all broadened (fa 2500 km s FWHM). 



lations of the structure of X-ray heated accretion disks indicate 
the existence of three sharply distinct layers near the surface: an 
X-ray heated, fully-ionized outer layer at the Compton tempera- 
ture (~ 10 7 — 10 8 K); a thin, highly-ionized layer just below that, 
with T ~ 10 6 K; and the usual, optically thick disk surface be- 
low that, with T ~ 10 4 -10 5 K (Nayakshin, Kazanas, & Kallman 
2000; Li, Gu, & Kahn 2001). Our measured temperatures and 
ionization states are thus consistent with those expected in the 
narrow, highly ionized layer just below the X-ray heated skin. 
For the observed O and Ne ion species to survive, the ioniza- 
tion parameter £ = L/nR 2 must be below 10 3 erg cm" 2 s -1 (Kall- 
man & McCray 1982), where n is the ion number density. This 
gives a rough lower limit of > 1 x 10 10 t/k pc cm for the radius 
of the line emitting region. This is consistent with an origin in 
the outer accretion disk, where X-ray heating effects dominate. 
We note that, for the inferred temperatures, plasma simulations 
with the XSTAR v2 code (Kallman & McCray 1982; Kallman 
& Bautista 2000) predict line emission from Si and S which we 
do not observe. This again points to an overabundance of Ne 
and O relative to solar values, as previously found by Angelini 
etal. (1995). 

4.2. Absorption Edges and the Nature of the Mass Donor 

Several of the observed photoelectric absorption edges are 
considerably stronger than expected for interstellar material, 
given the modest neutral hydrogen absorption determined from 



Lya measurements. We note that the most clearly detected 
strong edges are those of Ne and O, the same elements whose 
emission lines we discussed above. We thus conclude that the 
strong edges are due to absorption in cool, metal-rich material 
local to the source. (We may neglect local contributions to the 
hydrogen column density, as we know the donor in an ultracom- 
pact binary must be hydrogen-depleted; see Nelson et al. 1986.) 
Angelini et al. (1995) raised the possibility that 4U 1626-67 is 
embedded in the supernova remnant from the birth of the pul- 
sar, but as they point out, the absence of strong lines from Si, S, 
and Fe is inconsistent with this idea. Instead, we presume that 
this cool material originated in the accretion disk or the mass 
donor, so that its composition is then a clue to the nature of the 
companion. In the last column of Table 1, we calculate the col- 
umn density of material local to the source, assuming that the 
interstellar medium contributes A^h = 4.5 x 10 20 with standard 
elemental abundances. The local column density of the element 
with atomic number Z is then = TV™ -7V H A Z , where ' is 
the total measured column density of that element and Az is the 
number abundance of the element relative to hydrogen in the 
interstellar medium. 

We may compare our resulting abundance ratios with those 
expected for the possible companion types. The Ne/O abun- 
dance ratio of 0.22 ± 0.15 is much larger than expected for ei- 
ther a partially degenerate hydrogen-depleted star or a He white 
dwarf. On the other hand, 22 Ne is the next most abundant ele- 
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Table 4 
Emission Line Diagnostics 











T 


kT 


n e 


Ion 


H 


G 


R 


(10 6 K) 


(eV) 


(10 12 cm" 3 ) 


Ne IX (blue) 


4.85±0.73 


4.10±0.62 


0.33±0.08 


1.62±0.65 


140±56 


8.3±1.4 


Ne IX (red) 


7.87±0.91 


3.83±0.50 


0.29±0.07 


1.81±0.71 


156±61 


7.7±1.3 


O VII (blue) 


3.39±0.54 


2.68±0.40 


0.23±0.08 


0.98±0.35 


84±30 


0.62±0.09 


O VII (red) 


3.39±0.54 


2.18±0.35 


0.20±0.08 


0.91 ±0.65 


78±56 


0.85±0.12 



a // = h/r, G = (f +i)/r, and R = f ji, where r, i, and / are the fluxes of the resonance, intercom- 
bination, and forbidden components of the He-like triplet and h is the flux of the hydrogenic Lya 
line. 



ment in a C-O white dwarf, which has typical initial mass frac- 
tions of X c =X = 0.49 and X Ne = 0.02 (Segretain et al. 1994). 
This still gives a Ne/O ratio of only a few percent. However, 
a white dwarf that has cooled sufficiently will undergo crys- 
tallization (see, e.g., Shapiro & Teukolsky 1983), leading to 
chemical differentiation (Isern et al. 1991). In a C-O dwarf, 
this will cause all the Ne nuclei to settle into the inner few per- 
cent of the star and will also cause an O enhancement (Isern et 
al. 1991); detailed calculations predict X c « 0.19, X « 0.65, 
and ~ 0.16 (Segretain et al. 1994), consistent with our ob- 
served ratios. Note that since a mass of 0.02 M Q is required for 
a white dwarf secondary in this system, it is precisely the inner 
few percent of the star that is relevant, since the outer layers 
of the white dwarf would have been stripped away much ear- 
lier in the binary's accretion history. The C/O abundance ratio 
is poorly constrained by our observations, but the range of al- 
lowed values is consistent with the expected ratio in the core of 
C-O dwarf that has crystallized. 

If the feature near 9.5 A is in fact the Mg K-edge, then there 
is also a substantial overabundance of Mg relative to the so- 
lar value. The resulting Mg/O abundance ratio (0.17 ±0.13) 
would be quite high for a C-O dwarf, and instead suggests the 
possibility that the companion is an O-Ne-Mg dwarf. Although 
these stars have a typical initial composition of Xq = 0.72, 
X Ne = 0.25, and X Mg = 0.03 (see Gutierrez et al. 1996 and 
references therein), crystallization during cooling would likely 
lead to a substantially enhanced central concentration of Mg, 
by analogy to what is calculated for C-O dwarfs (Ogata et al. 
1993). Some models suggest that the central region might also 
contain some unburned 12 C from when the O-Ne-Mg dwarf is 
formed (Dominguez, Tornambe, & Isern 1994), which might 
account for the C overabundance suggested in our analysis. 

Taken together, the double-peaked emission lines and strong 
absorption edges of Ne and O in our HETGS spectrum thus 
point to the companion in 4U 1626-67 being the chemically 
fractionated core of a C-O or O-Ne-Mg crystallized white 
dwarf, the first example of such a donor among the low-mass 
X-ray binaries. In either case, the progenitors of both compo- 
nents of the binary must have been massive stars in an initially 
wide orbit in order to avoid a common envelope phase prior 
to the formation of a C-O or O-Ne-Mg degenerate core (see 



Levine et al. 1988 and references therein). The extremely low 
mass of the donor suggests an extended period of mass transfer. 
However, it is generally believed that prolonged mass transfer 
onto neutron stars leads to millisecond spin periods and weak 
(~ 10 8 G) magnetic field strengths (e.g., Bhattacharya & van 
den Heuvel 1991). By contrast, 4U 1626-67 has a slow (7.66 
s) spin period and a strong (~ 10 12 G) magnetic field. One 
possible explanation is that the neutron star is relatively young 
and that it formed from accretion-induced collapse of a mas- 
sive white dwarf (Taam & van den Heuvel 1986; Levine et al. 
1988; Verbunt et al. 1990). If that is the case, then the strong 
absorption edges we observe may actually arise from material 
expelled from this progenitor at the time of its collapse into a 
neutron star, rather than from the mass donor or the accretion 
disk. 

One further puzzle is to reconcile the observationally inferred 
mass transfer rate of M w 2 x 10" 10 M Q yr" 1 (Chakrabarty et al. 
1997; Chakrabarty 1998; Wang & Chakrabarty 2001) with the- 
oretical expectations. Mass transfer in a 42-min binary is driven 
by angular momentum losses via gravitational radiation (Ver- 
bunt & van den Heuvel 1995 and references therein). For a 0.02 
M donor, we expect only M « 3 x 1O~ U M yr" 1 . It is possi- 
ble that the currently observed rate does not reflect the long- 
term average rate. Indeed, we note that forM ~ 10~ n M Q yr" 1 , 
a C-O accretion disk may be subject to dwarf-nova-like ther- 
mal ionization instabilities (Menou, Perna, & Hernquist 2001). 
This may be true for O-Ne-Mg disks as well, although we are 
not aware of any explicit calculations. While 4U 1626-67 is 
generally considered a persistent source, its flux has been de- 
caying steadily for over a decade. It is possible that it is actually 
a very long-lived X-ray transient, and that its average M (cal- 
culated over many decades) is consistent with the theoretically 
predicted value. 
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